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Abstract 

This work analyzes the turbulent velocity distribution in proximity of the 
wall using the finite-scale Lyapunov theory just presented in previous works. 
This theory is here applied to the steady boundary layer under the hypothesis 
of moderate pressure gradient and fully developed flow along the streamwise 
direction. The analysis gives an equation for the velocities correlation and 
identifies the parameters of the expression of the average velocity through 
the statistical properties of the velocity correlation functions. In particular, 
the von Karman constant, theoretically calculated, is about 0.4, and the 
dimensionless Prandtl's length is in function of the Taylor-scale Reynolds 
number. The study provides the average velocity distribution and gives also 
the variation laws of the other variables, such as Taylor scale and Reynolds 
stress. The obtained results show that the finite-scale Lyapunov theory is 
adequate for studying the turbulence in the proximity of the wall. 
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1. Introduction 



The average velo city law near the wall is well known from many years 



von Karman J 193~ol) ) fo r what concerns the flow in smooth pipes f Nikuradse 



(119321 ) . iReichardtl ( 1l940l )) a nd in the cases of turbulent boundary layers with 



moder ate pressure gradient (iKlebanofil (Il952l ) , iKlebanofil (119551 ) , ISmith fc Walker 
(119581 ) ) . This law, usually expressed as 
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gives the average velocity in different regions near the wall, where y + = 
yUT/v and U + = U /Ut are the dimensionless normal coordinate and average 
velocity, and Ut = \/ v(dU/dy)o is the friction velocity. 

These expressions and the values of A, B, C and k, seem to be univer- 
sal properties of the flow, which do not depend on the Reynolds number. 
With reference to Fig. [TJ these equations, obtained through considerations 
of dimensional analysis and of self-similarity, hold under the hypothesis of 



fully d eveloped parallel flow along the streamwise direction x (jvon Karman 
(Il930h ). Specifically, Eq. JTJ, being the direct consequence of the wall bound- 



ary condition and of the definition of Ut, expresses the velocity distribution 
in the laminar sub-layer LL, a domain adjacent to the wall where the effects of 
the viscosity are dominant. Equation ([3]) holds in the turbulent region TR, a 
zone of non- homogeneous turbulence, where the Taylor scale Reynolds num- 
ber is high and var iable with y. I nto E q. (j3J), A; is the von Karman constant 
which according to lvon Karman! (119301 ) and in line with Eq. 03), can also be 
expressed as 



k = — lim 



d 2 U + 
dU + 



(4) 



where y£ defines the lower limit of the turbulent region. As far as Eq. (J2]) 
is concerned, this describes the velocity variations in the buffer layer BL, an 
intermediate zone between TR and LL in which viscous and inertia forces 
are comparable as order of magnitude. 

k and A are directly related to the velocity variations along y, whereas C 
and B express the order of magnitude of U with respect to the friction ve- 
locity in case of smooth wall. These constants ar e dimensionless free param- 
eters which can not be theoretically calculated (ILandau fc Lifshitzl (119591 )). 
ther efore several experiments d e aling with their de t ermin ation were carried 
out ( IFernholz fc Finleytl (119961 ). IZagarola fc Smitd (|1998l0. and s caling law 
similarities were proposed to justify Eq. (131) fearenblattl (jl993h ). These 
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Figure 1: Schematic of the boundary layer. 



constants can be also identified through the elaboration of t he resul t s of d i- 
rect numerical simulation s of the Navier-Stokes equations (jSpalartl ( 119881 ) . 
Fernholz &: Finleytl ( 119961 ) and references therein) . From the various sources 
of the literature, k = 0.36 4 0.44, C = 4.5 4 7.5, A ~ 5 4 6 and B ~ -5 4 -3. 

Another law of the wall, also valid under the conditions of parallel flow 
and self-similarity, and that seems to show properties of universality, pertains 
the Reynolds stress (u x u y ) 



(uv) + = g(y + ), (uv) + = (u x u y )/U 



(5) 



This behaves like g ~ y +3 very near the wall, and at about y + = 5, exhibits 
an inflection point whose coordinate provides the order of magnitude of the 
laminar sub-layer tickness. The value —(uv) + = 0.5, achieved in the buffer 
layer at about y + = 10 4 13, corresponds to the maximum turbulent energy 
produ c tion due to the mean flow, whereas far from the wall g ~ 1 (jHinze 
(ll975h . lTennekes fc Lumlevl <jl972h ). 

Although all these di strib utions seem to be universa l laws, according to 
Landau &: Lifshitzl (119591 ) and iFernholz &: Finleytl (119961 ) , there is not an ad- 
equate theory based on physical conjectures, other than the direct numerical 
simulations of the Navier Stokes equations, which leads to the calculation 
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of the free parameters of these laws. This represents the main motivation 
of the present work , whose pu r pose is to apply the finite scale Lyapunov 



theory, described in Ide Divitiisl (IDES 20111 ) . to determine the wall velocity 



distribution and its parameters. 



The firs t part of the work summarizes the main results of Ide Divitiis 



(lUES 20111 ). which deals with the steady homogeneous turbulence with uni- 
form velocity gradient, and then shows that such results can be extended 
also to the non-homogeneous turbulence with shear rate. We give a reason- 
able demonstration tha t the equation of the correlation function obtained by 



de Divitiisl (lUES 20111 ) can hold also in the case of non-homogeneous turbu- 
lence, in particular in the turbulent zone of the boundary layer. This result 
allows to express, in the turbulent region, the several quantities such as di- 
mensionless velocity gradient and turbulent kinetic energy in function of the 
local Taylor-scale Reynolds number, whereas in the laminar sub-layer and in 
the buffer region, adequate variation laws of the variables are assumed, which 
satisfy the wall boundary conditions and match the values of the variables 
in the adjacent domains. 

The analysis gives the von Karman constant in terms of the variables at 
the edge of the turbulent region, and provides a logarithmic velocity profile 
in TR, different from Eq. ([3]), which exhibits free parameters (k included). 
The knowledge of the statistical properties of the spanwise correlation func- 
tion leads to the estimation of the variables at y+ and thus to identify k 
and the other parameters. Therefore, the von Karman is here theoretically 
obtained, resulting to be about 0.4, in line with experiments and numerical 
simulations, and the other results, in agreement with the several data from 
the literature, show that the finite scale Lyapunov theory can be an adequate 
tool for studying the wall turbulence. 



2. Resume 



This section summarizes the main results of lde Divitiisl (IDES 20 111 ) which 
regards the steady homogeneous turbulence in the presence of an average ve- 
locity g radient V X U. Th ere, the author, applying the Finite-scale Lyapunov 
theory (Ide Divitiisl ( 120101 )) and the Liouville theorem, proposes the following 
evolution equation of the pair distribution function of fluid velocities in 
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case of arbitrary flow 



dp® 



+ Vx F( 2 ) • v + V X ,F( 2 ) • v' = A(r) (f o (2) - F^) - J D 



(6) 



(2) 

where F is the pair distribution function of the isotropic turbulence which 
exhibits the same momentum and kinetic energy of F^ 2 \ whereas —Jd repre- 
sents the rate of caused by the rate of the turbulent kinetic energy. A(r) is 
the finite-scale Lyapunov exponent associated to the finite-scale r = |x' — x|, 
and v and v' are the fluid velociti es calculated at x and x' respectively. 
In case of homogeneous turbulence, Ide Divitiisl ( jlJES 20111 1 shows that the 
steady distribution function reasonably tends to a quantity which depends 

(2) 

upon Fq , V X U and A 



d4 2) dU Sv + dF, 

dv', dx 



,(2) 



A(r) I dvj dx 



m 



/' 



The velocity correlation tensor R ki = (ujw') is then calculated, by definition 



Rki 



F^u k u' d 3 u d 3 u' 



p U dUk n + dUi R 



(8) 



where R^ is the sec ond order velocity correl a tion t e nsor assoc i ated t o the 
isotropic turbulence ( von Karman fc Howarth ( 1938 ). Batchelor f 19531 )) 



Rkio(r) 



u 



(f-g) 



+ 



>ki 



(9) 



being / and g — f + 1/2 r df /dr longitudinal and lateral velocity correlation 
functions, respectively, and u = v — U is the fluctuating velocity. For r = 0, 
Eq. (jBJ) provides the expression of the Reynolds stresses in function of V X U 



{u k Ui) 



u 



1 (dU k dUi 

Ohi I — h 



A V dxi dxf. 



(10) 



Equation f JTQj) is a Boussinesq closure of the Reynolds stress, being u = 
(uiUi)/3, (.) indicates the average calculated on the ensamble of the fluid 
velocity, and A = A (0) is the maxima l Lyapunov exponent. The Lyapunov 
theory presented in Ide Divitiisl ( 120101 ) shows that A = u/Xt, being At the 
Taylor scale. 
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Next, the conditi on of steady flow leads to the following ordinary differ- 
ential equation for / Ide Divitiia ( llJES 201 ll ) 



1 - / df 2 fd 2 f 4df\ 10 fr 



2 rff R T \ df 2 fdf J R T y/2(l - f) 

whose boundary conditions can be reduced to the following conditions of / 
in the origin r = 

/(0) = 1, ^ = (12) 

where f = r/Ay, d 2 f/df 2 (0) = —1, and Rt = u\t/v is the Taylor scale 
Reynolds number. Thus, Eqs. (ITT|) and ( IT2|) express an initial condition 
problem whose initial condition is given by Eqs. ( !T2l . The same steady con- 
dition gives the relationship between V X U, A and Rt, which represents the 
equation of the kinetic energy in the case of steady homogeneous turbulence 



S_ 15_ 

A 2 ~ R T 



(13) 



where 5* is related to V X U 



dxk V ® x % ® x k 

Th e tensor R i k is then cal culated with Eq. §5§ and ([2]), by means of /. 



de Divitiisl ( llJES 201 ll ) studies some of the properties of Eqs. fTTT]) - ( p~2|) . 
and shows that their solutions behave like / — 1 ~ r in a given interval of r. 
Accordingly, the statistical moments of velocity difference are (Au n ) ~ r n l 2 
for n < 5, and the energy spectrum E(k) w fc -2 in the inertial subrange. 

Equation (jH]) was obtained for an arbitrary flow, whereas the other equa- 
tions were derived under the assumption of steady homogeneous turbulence 
with a given average velocity gradient. The successive sections show that 
these equations can be applied also in the case of non-homogeneous turbu- 
lence where At and u vary with the space coordinates. Thus, the method is 
here applied to the turbulent region of the fully developed boundary layer, 
where the several quantities vary with the wall normal coordinate. 
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3. Analysis 



An evolution equation for the velocity correlation is determined, in the 
case of non-homogeneous turbulence with a nonzero average velocity gradi- 
ent. 

The fluid velocity, measured in the reference frame 5i, is v = U + u, where 



U = (U x , U y , U z ) and u 



u z ) are, average and fluctuating velocity, 



respectively. The velocity correlation tensor is R^ = (uiUj), where Ui and u'j 
are the velocity components of u calculated at x and x' = x + r, and r is 
the separation distance. As the analysis is finalized to the description of the 
turbulent region of the developed boundary layer, and since there the effects 
of the spatial variation s of V X U are orders of magnitude muc h smaller than 
those caused by V X U f von Karman ( 1930l ). Nikuradse ( 19321 )). the velocity 
gradient is assumed to be a function of x alone, being V X U = V X U'. 

In order to determine the evolution equation of R^, the Navier-Stokes 
equations are written for the fluctuating velocity in the points x and x'. 
The evolution equation of Rij is determined by multiplying first and second 
equation by u'j and u i: respectively, summing the so obtained equations, and 
calcul a ting the ay e rage o n the statistical ensemble (jvon Karman fc Howarth 
fll938h . lBatchelori dl953h ) 



dR 



dt 
being 



, 2 dllj dRij 

Mj + Pij + 2z/V R^ — -^-Rkj — -^—Rik + a _ [Uk 



dx k 



dx k 



dr k 



U' k ) + T l3 (lh) 



T i:j {x, r) 



_d_ 



[uiu'^Uk - u' k )) , 



1 d(pu' j ) 
P dri 



(16) 



and p is the fluctuating pressure. The quantities of Eq. f|T5l) . which in turn 
depend on r, due to non-homogeneity depend also on x, and Eq. ( |T5i) incudes 
an additiona l term with resp ect to the homogeneous turbulence, represented 
by r -(x, r) flOberlack! (Il997h ) 



r^(x,r) 



_ d Rij 
-2v- 3 



1 d(pu'j) 
p dxi 



_d_ 

dxi 



(UiUkU'j) 



1 d(p'ui) 



p drj 



+ v 



d 2 Rjj 
dxi-dxi 



(17) 



dr k dx k 

which provides the non-homogeneity of the different terms of correlation. 
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Making the trace of Eq. (IT5|) . we obtain the following scalar equation 
dR 



being Rf k is the symmetric part of Rik and 



R 



Ri 



r 



(18) 



(19) 



i?(x, 0) gives the turbu lent kinetic energy, H = Tg p r ovides the mecha- 
nism of energy cascade ( von Karman fc Howarthl ( 1938T ) . iBatchelorl ( 19531 )). 
Pn = expresses the fluid incompressibility whereas T arises from the non- 
homogeneity of the flow. 

Now, a scalar equation for describing the main properties of the velocity 
correlation is determined. To this end, Rij, H and U are decomposed into 
an even function of r = |r| (here called spherical part), plus the remaining 
term: 

Rij(x, r) = Rij(-x, r) + ARij(x, r) 



#(x, r) = H (x, r) + A#(x, r) 



(20) 



U'-U = U(x,r) + AU(x, r) 
where F is the spherical part of the generic quantity F, defined as 

F(x, r) = - (F(x, r, 0, 0) + F(x, 0, r, 0) + F(x, 0, 0, r)) 
o 



+ (F(x, -r, 0, 0) + F(x, 0, -r, 0) + F(x, 0, 0, -r)) 
o 



(21) 



and ARij(x, 0)=Aif(x, 0) = 0. Therefore, the Fourier transform of R identi- 
fies the part of the energy spectrum depending upon k 2 . As the consequence 
of this decomposition, R and AR satisfy the equations 



dR _H f d 2 R 2dR 
dt 2 \ dr 2 r dr 



(22) 



(23) 



in which Eq. (123 j) is obtained as the difference between Eqs. (1181) and (I22p . 
and 



G 



m 

dx k 



R s 



kl 



Go, AG 



ki 



. ( dR . 



U'k) - Ar (24) 



and G represents the spherical part of —dR/dr k {Uk — U' k ) — Y. It is worth 
to remark th a t Eg. ([221 formally coincides with the equation obtained 
de Divitiisl ( IDES 20111 ). with the difference that here, because of non- 
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homogeneity, the quantities appearing into Eq. (|22|) depend also on x. 

The turbulence is here studied using Eq. ( 122|) alone, whereas i?y(x, r) will 
be determined in function of V X U, by means of a proper statistical analysis 
of the two-points velocity correlation. 



4. Pair distribution function 

This section analyses the non-homogeneous turbulence through the pair 
distribution function, taking into account that At and (uiUj) vary with the 
spatial coordinates, whereas V X U is considered to be an assigned quantity. 
To study this, consider now the pair distribution function of the fluid velocity 

( v, v'; x, x') = Ff ( v, V; x, x') + ^ (v, v'; x, x') (25) 

where obeys to Eq. ()6]), and (j^ 2 \ representing the deviation from 
the isotropi c turbulenc e, satisfies, at each instant, the following equations 
fide Divitiisl fllJES 201lh ) 



{2) du 3 du' 3 = 0, 
j j <f) {2) u du 3 du' 3 = 0, (26) 

J V J v' 



i (2) u- u du 3 du' 3 = 0. 



Equations (126]) state that momentum and kinetic energy associated to F 
and F^ are equal each other, respectively. The analysis supposes also that 
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all the dimensionless statistical moments of F are constant in all the points 



of the fluid domain, therefore the functional form of _F is assumed to be 



F (2) (v,v';x,x') 



,(2) 



Ufx) v' 



,(2) 


Ufx') 



U X 



u(x') 



(27) 



and its gradients are calculated in functions of the spatial derivatives of U 
and u 



dF, 



(2) 



dF, 



(2) 



dxi 



dvi V dx 



dUj Uj du 



u dxh 



dF, 



(2) 



<9F (2) fdUj u'j du' 



dxi 



dv'i V dx 



+ 



u' dx'i 



T (28) 



where, as before, V X U = V X 'U. Substituting Eqs.fl25l) and ( 1281 into Eq. 
(E]), we obtain the following relationship between F^ and 



A0 



(2) 



J 



D 



dFf 

at 



+ 



2) 



dt 



dx p 



v v + 



dx' p !i 



+ 



dF^ v p Uj du 
dvj u dxk 



+ 



dF (2) v'pu'j du' 
dv'j u' dx' k 



(29) 



+ 



dF, 



(2) 



<9f ,■ 



9F f 



(2) 



M Vp 



dx n 



As the last term of Eq. f )29|) identically satisfies Eqs. ff26|) . these latter are 
written in the form 



1 
u 
u • u 

^(2) „'„' 



9iT(2) ^ ^(2)^ | ^(2)^ 

<9t> 7 ' It tfofc 



<9t 



+ 



dx n Vp 



dx' v v 



(30) 



dF^> v'pu'j du' 
dv'j u' dx' k 



J D ) du 6 du' 6 = 



A sufficient condition for satisfying these five equations is that the integrand 
of Eq. ( 1301) vanishes 



&p(2) ^(2) d( p(2) f dF w VpU . du dF m v > pU >. du , 

+ — — v v + ^ ■ v„ — — —TT~r + Jd 



dt 



dx p "'' 



dx' p ->' 



dvj u dxk 
10 



dv'j u' dx' k 



0(31) 



Assuming that Eq. (13~Tj) is true, and taking into account Eq. (|2T?j) . -F^ 2 ^ is 
given by Eq. (j7j) which does not depend on the gradient of u and w', as in 
the case of the homogeneous turbulence. Really, changes starting from 
an arbitrary initial condition, therefore Eq. (J7J) represents an approxima- 
tion which can be considered to be valid far from the initial condition. As 
the consequence, also the spherical part of the correlation tensor, R\j, here 
obtained 



R 



ki 



R 



ik 



- (V + —r 

3 \ dr 



ki 



dU k | dUj 
dxi dxk 



(32) 



coincides with the expression given in Ide Divitiia ( IDES 20111 ). Thus, sub- 
sti tuting Eq. (132]) and int o Eq. ( 1221) and following the analytical procedure 
of Ide Divitiisl ( IIJES 20111 ). we found that / obeys to Eq. ffTTT) also in the 
present case, and A is related to VxU and Rj through Eq. ({TBI . The dif- 
ference with respect to Ide Divitiisl (IIJES 20111 ) is that, here the turbulence 
is non-homogeneous, thus Rt = Rt(x), and / also depends on x, being 
/ = /(x,r). 



5. Analysis of the boundary layer 

Here, the steady turbulent boundary layer with a moderate pressure gra- 
dient is analyzed, assuming that the flow is fully developed along the stream- 
wise direction. To this purpose, return to Fig. [IJ and consider only the de- 
veloped region of the flow. In the figure, is the wall frame of reference, x 
and y are, respectively, the streamwise direction and the coordinate normal 
to the wall, whereas z is the spanwise coordinate. 

In the developed region, u, (u x u y ) and At change with y, and in the 
laminar sublayer U x and u are both about proportional to y. In LL, BL and 
at the beginning of TR, the correlation scale of velocity is proportional to 
the distance from the wall, and vanishes for y — 0, being 

*-($)." + - (33) 

where (d\T/dy) = 0(1), whereas the velocity fluctuations follow the Navier- 
Stokes equations and satisfy the wall boundary conditions 

u x = d -^(0)y + ..., % = ^(0)l/ 2 + ..., u z = ^(0)y + ... (34) 
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In case of fully developed flow along x, that is, parallel flow assumption 
(d/dy > >> d/dx), t he co ntinuity and momentum equations of the mean 
flow are flSchlichtind (12004 )) 



9Uy , , 

■ai = °- < 35) 

dy x v p dx dy 2 

(36) 

9 2 i ap () 

dy y p dy 

where P is the average pressure, whereas the equation of the turbulent kinetic 
energy reads as 

d I ( p ujUj \ \ , > dU x I duj du-, 



^^ + ^r)) + ^W +, '\^) = (37) 

The non-homogeneity is responsible for the first term of Eq. (!37j) . whereas 
second and third terms represent, respectively, the energy production due to 
the average motion, and the dissipation. 

Because of the parallel flow assumption, in all these equations U x , (uiUj) 
are functions of y alone. From Eq. (135]) and taking into account the boundary 
conditions, U y = 0, whereas Eqs. (13"61) give the following first integrals 

E^l = Fx + H- (u 2 y ) (38) 



(u x u y ) - = -Fy - If? (39) 

being F = 1/p dP/dx and if is a proper constant proportional to the average 
pressure at x — 0. Into Eq. ( 1391) . the boundary layer approximation and the 
hypothesis of moderate average pressure gradient along x provide that 11% » 
\Fy \ in all the regions, therefore introducing the dimensionless variables U + = 
U x /Ut and y + = yUx/v, Eq. ( l39|) reads as 

^ = (40) 
dy+ 1 + R T 
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Equation f l4"U|) is assumed to describe the average flow with moderate pressure 
gradient in the three regions of the boundary layers, when the parallel flow 
hypothesis is verified. 



5.1. The Turbulent Region 

This section studies the distribution of the different variables in the tur- 
bulent region TR, by means of the analysis seen in the sections [3] and HJ 

First, observe that, in case of fully developed parallel flow, the effects 
of non-homogeneity in TR are much smaller than tho se related to energy 



production and dissipation (ITennekes fc Lumleyl ( 119721 )). thus the first term 



of Eq. fl37|) is here neglected with respect to the other ones, and Eq. f fT3|) 
is recovered. This approximation, in agreement with the analysis seen in 
sect. HI states that the kinetic energy production is balanced only by the 
dissipation in TR, and allows to express the several variables in terms of the 
local value of Rt (or dU x /dy). 

In order to calculate u + = u/Ut, dU x /dy is eliminated between Eqs. ( 139|) 
and ([T3D, where S = (dU x /dy) 2 

o3/4 

u + = K T (41) 

15VV1 + R T 
Being Rt = u + X^, also A J is in terms of Rt 

A+ = (lhR T ) l/ yi + R T (42) 

As these equations arise from Eq. (1131) which holds in the turbulent region, 
Eqs. (l4"Tj) and (14"2"1) describe the variations of u and At only in TR. 

From the comparison between Eqs. f|4T)l) and (@J, (dRT/dy + ) e identifies 
the von Karman constant, here expressed taking into account that Rt = 

where the subscript e denotes the values calculated at the edge of TR. Substi- 
tuting Eqs. (HT!) and (H2|) into Eq. (1431) . we obtain the von Karman constant 
in function of the variables at y+ 



dX T '\ a/1 + i?Te D 3/4 
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It is worth to remark that this estimation of k requires the knowledge of 
Rxe and of (d\ T ~/dy + ) e , whereas does not need the assumption that U + is 
represented by a logarithmic profile. 

Next, the dimensionless Prandtl's mixing scale is calculated from the 
definition of Prandtl's mixing length l p 



- {U X Uy} 

and Eqs. (110]) and (fB]). This is 



dy 



W l * (45) 



i; = ^/R T (l + R T ) (46) 

being Z+ ~ R T for R T » 1, and its derivative calculated at yf is also 
expressed in terms of Rx e 



dy+J e 2^/R Tt (l + R Te ) 

In view of Eq. f j44l) . this derivative can be also expressed in function of 

dX J / dy 

d£\ _ l + 2R Te f R Te \ 1/A ( d\+\ 
dy+) e l + 3i? Te V 15 J \dy + J e 1 1 

This expression gives the link between the variations of the Taylor scale and 
of the Prandtl's length at the border of the turbulent region. 

Following Eqs. ( 14 ip and (I42p . u + and A J are functions of y + through 
the local value of Rt (or of dU + /dy + ), therefore, the distribution of such 
quantities along y + require the knowledge of the function Rt = Rt{i) + )- 
This latter can be expressed as 

R T (y + ) = R Te + (y + - yt) + 0{y+- ytf (49) 



UOl'O « ) ( // ' — /// 1 

other terms, dll + /dy + is 
dU + 



Now, in a range of y + where O (y + — y£) is negligible with respect to the 

(50) 



dv+ ^ R -<W)y-^ 
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U + exhibits there logarithmic law, obtained integrating Eq. ( 15 Op from y£ to 

y + 

^M l+RT \ + + Z- Vt) )^ 

being = U + (y+). This law is defined as soon as the parameters y+, k 
and are known. Equation (I5TT) differs from the classical expression (J3j 
and formally tends to Eq. ([3]) when y + — > oo. Therefore, Eq. (l5"Tj) can give 
values of t/ + sizably different from ([3]) for small y + . The comparison between 
these equations, for y + — > oo identifies C in terms of the variables at y e 

C = U+ + Un(-±-) (52) 



fc V 1 + R Te 



As far as the Reynolds stress is concerned, it is expressed in function of 
R T through Eqs. (EHD and (SO]) 

<™> + -^-ifk < 53 > 

Observe that, the wall boundary conditions state that (uv) + y +3 near 
the wall, whereas Eq. (|53|) gives (uv) + m R T « y +2 . This disagreement 
is due to the fact that the expression of (uv) + has not been derived from 
the correlation equation with r = 0, but arises from Eq. (|39l) eliminating 
dll x /dy. This implies that Eq. fl53|) holds only in TR and BL, whereas in 
the laminar sub-layer a proper matching condition must be applied. 



5.2. Matching Turbulent region - buffer layer 

With reference to Fig. [TJ the domains LL and BL constitute SL, a zone 
between wall and turbulent region. There, due to the presence of the wall, 
the analysis of sections H] and [3J can not be applied. Therefore, the mean 
variables in SL are expressed in function of y + , taking into account the 
boundary conditions and that At follows Eq. (133|) . where it is assumed 
(d\^/dy + ) e = (d\^/dy + )o. As the consequence, Rt, u + and are supposed 
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to vary in SL according to 

Rt = w^Aj, 

_2 



, y + + C u y^ . 
" = u: — — — , (0 < y^ < y e 



vt + C u yt ^> 
A? = A+ e ^, (0 <!/+<?+), 

tie 

where C u is a constant given by 

C. = (55) 

Equations (1541) and ( 1551) provide the matching condition between SL and TR. 
Specifically, Eqs. (I5l| state that u + , Rt and are continuous functions for 
y + = vt i whereas Eq. fl55l) gives there the continuity of their derivatives. The 
average velocity is then calculated by quadrature, substituting the expression 
of Rt given by Eqs. fl54|) . into Eq. (j40p . integrating this latter from to y + , 
with C/+(0) = 

^^riTW)^ <56) 

where U + (y + ) ~ y + for small y + . The matching between SL and TR gives 
the value of . 

^ - f + TTW)^ (57) 

Thus, Eq. ( 1561) and ( ISTi) establish that both £7 + and dU + /dy + are continuous 
for y + = y+ 



5.3. Matching buffer layer - laminar sub-layer 

As previously seen, the expression ( |53l) of the Reynolds stress, valid in 
BL and TR, can not be applied in the laminar region. Since (uv) + ~ y +3 
near the wall, the Reynolds stress in LL is approximated by 

(uv)+ = (uv)t y ll + ° uvy+ + l (0<y + <yt) (58) 
y* + ^uvV* 
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being C uv a constant 

c ^ = 1 3-(d\n(uv) + /dy + )*yt 
yt (dln(uv) + /dy + )*yt - 4 

Equations f )58|) and fl59|) establish that the Reynolds stress and its deriva- 
tive are continuous for y + = where the subscript * indicates the value 
calculated at y+. This latter, obtained as the inflection point of (uv) + (i.e. 
d 2 (uv) + /dy +2 = 0) according to Eqs. flo^]) and (Ell), gives the dimensionless 
tickness of the laminar sub-layer. 



6. Identification of the velocity law free parameters 

The definition of the velocity law, requires the knowledge of the parame- 
ters which appear into Eqs. (jUJ) and floTj) . In particular, the determination 
of k needs the values of Rx e and (d\^,/dy + ) e . To identify these latter, we 
will proceed as follow. 

First, observe that (uiUj) is a symmetric tensor which can be obtained 
from the diagonal tensor 



(u u ) 











(vl) 











(60) 



through an opportune rotation around y = yo, being (-Uq), (vq) and (w^) the 
velocity components standard deviations in the canonical frame, and $ is the 
angle of this rotation. Therefore, u and (u x u y ) are in terms of the elements 
of (u u ) 



u 2 = M) + (v 2 ) + K 2 » 



(U X Uy) = ((U 2 Q ) ~ (V%)) 



sin 2d 



(61) 



Because of the parallel flow assumption and taking into account Eqs. ([10 
and ( l6il . u and (u x u y ) are related each other in such a way that 



U > \ {U X Uy) 



(62) 



17 



Hence, Eq. (I13p implies that 



dU x , 

Rt > 15 or — — \t > u 
dy 



(63) 



in the turbulent region. This limitation identifies the minimum value of Rt 
in TR, which is assumed to be 



R 



Te 



15. 



(64) 



To estimate (d\ T ~/dy + ) e , consider first the spanwise correlation function of 
the streamwise velocity components (that is Rn(r z ) = (u x (x, y, z)u x (x, y,z + 
r z ))). This can be calculated with Eq. (jSJ), once / is known through Eqs. 
(TTT1) . Since dU/dy leads to the development of coherent structures in the 
fluid, similar to streaks and ca used by the stre tching of the vortex lines 
along the streamwise direction f Kim et all (jl990h ). we expect that Ru(r z ) 
intersects the horizontal axis and remains negative for r z — > oo. This im- 
plies a wide distribution of spacings between the differ ent streaky structures, 
whose mean value depends on Rt (IKim et al.l ( 119901 )) . Accordingly, there 
exists a minimum spanwise distance r* such that, for r z 6 [0, r*], R\i{r z ) 
gives the ne cessary i nform ations to describe the stati stical pr o pertie s of the 
correlation ( Ventsel ( 1973 )). According to the theory ( Ventsell ( 1973 )). if Rn 
monotonically tends to zero as r z — > oo and Ru < 0, r* can be estimated as 
the distance at which Rn is negative and exhibits an inflection point (second 
inflection point of the curve, see Fig. [2]) 



R n (rl) < 0, 



d 2 R 



ii 



(r*) = 0, Condition 1 



(65) 



Alternatively, r* can be estimated as the intersection between the osculating 
parabola in the point f z where Ru = min, and the horizontal axis, that is 



Rn(r z ) + - 



ld 2 R 



2 dr 



r z ){r z 



0, Condition 2 



(66) 



Since Rn(r z ) is expressed in function of / through Eq. (jSJ), this is related 
to the correlation functions associated to the other directions, thus r* is 
representative also for the other coordinates x and y. As the result, the 
distance from the wall of a point of TR must be always greater than r* (i.e. 
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r* < y). Hence, it is reasonable to assume that, at the edge of the turbulent 
domain 

Ve = K 

(67) 

R Te = 15 



In order to calculate Rn, f is first obtained by solving Eqs. flTTl)-ffT2l) 
which corre spond to the following Cauchy's initial condition problem (jde Divitiisl 
f llJES 201l( l ) 



df 



dF 5/r 




(68) 



/(0) = 1, F(0) = 
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Hence, R\i{r z ) is obtained with Eq. (jSJ), and r* is calculated for both the 
conditions (1B"5|) and ( 16"6"j) . As the result, y+ is given by 

Jk- = J!*, where i? T = i? Te = 15 (69) 
Are Are 

and (d\^/dy + ) e is determined according to Eqs. (15"4j) 

=*| 



7. Results and Discussion 



As the calculation of and dX^/dy + needs the knowledge of the sta- 
tistical properties of the velocity correlation, i2 n is calculated for different 
values of Rt- To this purpose, / was first determined by solving numerically 
Eqs. (|68p . by means of the fourth-order Runge-Kutta method. The calcula- 
tion was carried out for Rt = 15, 20, 40, 60 and 80, where these Reynolds 
numbers correspond to several distances from the wall in the turbulent re- 
gion. To obtain a good accuracy of the solutions, the step of integration 
is chosen to be equal t o 1/40 of the est imated Kolmogorov scale Ik, where 
Ik/X t = l/15V 4 /y/Ifa flBatchelori dl953h l The results are given in Fig. E2(a) 
which shows / in terms of r, where the bold line represents the correlation 
function at the edge of TR (Rj< = 15). This latter exhibits the ratio (integral 
scale) / (Taylor scale) quite similar to that of a gaussian centered in the origin, 
whereas in the other cases , this ratio increases with Rt, in agreement with 
the analysis of Ide Divitiia (llJES 20111 ). A more detailed analysis about the 
changing o f f and of t he correspo nding energy spectrum E(n) with Rt is 
reported in Ide Divitiia ( IIJES 201ll ). 

The spanwise correlation function Ru(r z ) is then calculated with Eq. OH]), 
and is represented in Fig. [3] (b) for the same values of Rt- From these data, 
the edge of TR, y+ = r*, is calculated for both the conditions 1 and 2. 

The von Karman constant and all the others parameters are shown in 
the table [TJ The table reports also the values of the free parameters A, B 
and C associated to the classical wall laws flTJ-flHJ), which are here calculated 
through the identification with Eqs. ( 151 j) and ( I56p . In particular, C is cal- 
culated with Eq. ( 152|) . whereas A, which pertains the velocity law in the 
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buffer layer, is identified as the slope dU + (y + )/d\ny + of Eq. (1551) where 
d 2 U + (y + ) /d\ny +2 = 0, and B is consequentely determined. The values of k 
and C are compared with the results given by different authors. 

For both the conditions, y+ ~ 50 and U£ ~ 15 and this c orresp onds 
to a difference w ith respect to the classical data of iNikuradsd ( 119321 ) and 
Reichardtl (jl940h . which is less than 2 %. As far as the Prandtl's length is 
concerned, according to Eq. (HEJ), it varies quite similarly to Rt, and its 
derivative, calculated for y + = y+, is slightly greater than k (see Eq. ( [47]) ) 



dl+ 



1.34782... 



dX 



r 



dy" 



1.00052 k 



(71) 



Its specific values, shown in the table, are in excellent agreement with the 
classical results. For what concerns A, B and C , their values agree quite well 
the experiments, expecially for what concerns C and A. 

Once the free parameters are identified, Rt is calculated in function of 
y + through Eq. 012]) with 0(y + - y+) 2 ->■ and (154]) . for both the con- 
ditions (see Fig. H]), and this corresponds to consider only the logarithmic 
profile of U + . Figure [5] shows U + (y + ) calculated with Eq. (15T|) and (156|) . 
It is apparent that the results agree v ery well w i th the experimental dat a 
of different authors (iNikuradsel (Il932h. iKlebanoffj ( Il952h. iReichardtl jl940h l 
and with the formulas proposed by SpaldingT ( 1961 ) and Musker ( 1979 ). with 
an error which does not never exceed 5 % in the interval of y + of the figure. 
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In particular, the condition 1 (Eq. ( 1651) ) provides a maximum difference be- 
tween present results and the data of Nikuradse and Reichardt, less than 3 
% in this interval. 

The other variables are represented in Fig. [6j The dimensionless Taylor 
scale (Fig. [6] a), linear in the laminar sub-layer and in the buffer region, 
remains a rising function of y + in the turbulent zone and exhibits, there, a 
slope which decreases slightly with y + , whereas the Lyapunov exponent (Fig. 
|6]b), defined only in TR, is represented by a monotonically decreasing func- 
tion of y + . This latter, being proportional to the square root of the turbulent 
dissipation rate (A = u/Jvr), agrees with the data of the different experiments 
( iFernholz fc Finleytl ( 119961 )). at least where U + exhibits logarithmic profile. 

The square root of the kinetic energy and the Reynolds stress are shown 
in Fig. |6] (c) and (d). These vary in TR following Eqs. f|4T]) - fl53|) and are 
monotonically rising functions of y + . As seen, u ~ y + and (uv) + w y +3 in the 
laminar sub-layer, whereas the Reynolds stress calculated with Eq. f )53|) in 
the buffer region, shows an inflection point y~£ which represents the separation 
element between LL and BL. For y + > y~£, — (uv) + rises with y + until to reach 
the turbulent region where is about constant and equal to the unity. More 
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Figure 5: Dimensionless average velocity profile: SL dashed lines, TR continuous lines 



in detail, Fig. [7] shows —{uv) + in an enlarged region which includes LL and 
part of BL. The distance represented by F\ and F 2 is y+ = 5.734 and 
6.097, in line with the order of magnitude of the laminar sub-layer tickness, 
and this is achieved at about — {uv) + = 0.24 in the two cases. Next, the value 
— (uv) + = 0.5 is obtained in the buffer layer for y + ~ 10.5 and 11.2, v a lues i n 
very good ag reement with the classical results (ITennekes fc Lumleyi (119721 ). 



Hinzd (119751 ) ). This last condition corresponds also to the maximum of the 
kinetic energy production in the buffer layer, as shown in Fig. [HJ In this 
figure, the variations of the rate of turbulent kinetic energy production due 
to the mean flow a re shown in LL and BL. The shape of the d iagrams follows 
the classical data ( Tennekes k, Lumley ( 19721 ). Hinze ( 1975 ) and references 
therein), and the values of (uv) + , (uv) + dU + 1 dy + and of the corresponding 
y + , are in excelle nt agreement with the data of ITennekes fe Lumleyi (Il972l ) 
andlHmzd ril97fih . 

It is worth to remark that, although the monotonic trend of u and (uv) + 
can contrast some experiments which can give non-monotonic variations of 
these variables, the values of such quantities and the corresponding y + , 
are co mparable with those of the several experiments (IFernholz fc Finleyt 
(119961 ) ) . This discrepancy can be due the fact that, here, the effects of 
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Figure 6: Distribution of the dimensionless variables in the boundary layer: SL dashed 
lines, TR continuous lines, (a) Taylor scale, (b) Maximal Lyapunov exponent, (c) r.m.s. 
of fluctuating velocity, (d) Reynolds stress. 



non-homogeneity of (u y (p/p + UjUj/2)) on the kinetic energy equation are 
neglected in TR, and that the average velocity is analyzed only in the loga- 
rithmic range. 



8. Conclusions 

This work analyzes the turbulent wall laws through the Lyapunov theory 
of finite scale. The results, valid for fully developed flow with moderate pres- 
sure gradient, are subjected to the hypothesis that in the turbulent region, 
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Figure 7: Dimensionless Reynolds stress in the laminar sub-layer and in the buffer region 



the energy production due to the average flow is balanced only by the dissi- 
pation rate, whereas the non-homogeneity of (u y (p/ p + UjUj/2)) is neglected. 

The free parameters of the velocity law, here theoretically calculated 
through the statistical properties of the velocity correlation functions, and 
the wall laws, are in very good agreement with the literature. In particular: 

• The von Karman constant, theoretically identified as k = (dRx / dy + ) e ~ 
0.4 -=-0.42, depends on the scale of the spanwise velocity correlation and 
does not requires the assumption of the logarithmic velocity profile. 

• The average velocity law and the distributions of the other dimen- 
sionless quantities such as kinetic energy and Reynolds stress in the 
boundary layer, agree -or at least are comparable- with experiments 
and direct simulations. 

These results, which represen t a further application of the analysis pre- 



sented in Ide Divitiisl ( 1IJES 20111 ) , show that the finite scale Lyapunov anal- 



ysis can be an adequate theory to explain the wall turbulence. 
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Table 1: Velocity law parameters and comparison of the results. P.R. as for "Present 
Result" . 

Parameter P. R. P. R. 

Condition 1 Condition 2 Nikuradse Zagarola Barenblatt Smith Fernholz 
Eq. gg Eq. (gg) 

k 0.4233 0.3984 0.4 0.41 0.425 0.461 0.4 

C 6.2445 6.4822 5.5 5.2 6.79 7.13 5.1 

A 5.2384 5.5652 

B -4.1226 -4.7166 



r* 

z 

dX 



dy- 



3.18233 3.38085 



£ 0.31423 0.29578 



dl+ , 

2r 0.4235 0.3986 



dy^ 

yt 49.3005 52.3760 

0.96825 0.96825 

-(uv)+ 0.93750 0.93750 



U+ 14.8250 15.7498 
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